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Holographic characterization of imperfect colloidal spheres 
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We demonstrate precise measurements of the size and refractive index of individual dimpled colloidal spheres 
using holographic characterization techniques developed for ideal spheres. 


Holographic snapshots of colloidal spheres can be in¬ 
terpreted with the Lorenz-Mie theory of light scattering 
to measure an individual sphere’s three-dimensional posi¬ 
tion, size and refractive index 1 . When applied to dielec¬ 
tric spheres with near-ideal sphericity and smoothness, 
this technique yields nanometer-scale precision for the 
position and radius 1 4 and part-per-thousand precision 
for the refractive index 1,5 . Here, we demonstrate that 
this technique yields similarly precise and meaningful re¬ 
sults for imperfect spheres, provided that their deviation 
from sphericity is not too pronounced. 

Our model system for this work consists of monodis- 
perse colloidal particles synthesized through emulsion 
polymerization of 3-methacryloxypropyl trimethoxysi- 
lane (TPM) 6 . Depending on how they are made, these 
particles can take the form of spheres, as shown in 
Fig. 1(a), or dimpled spheres, as shown in Fig. 1(b), in 
which both the sphere radius and dimple dimensions are 
drawn from narrow distributions. 

The two types of particles are synthesized through sim¬ 
ilar pathways. TPM oil, which ordinarily is insoluble in 
water, undergoes hydrolysis in a basic environment (pH 
> 9) and becomes water soluble. Solubilized monomers 
then form insoluble oligomers, which condense into spher¬ 
ical droplets. Once the droplets are fully grown, they are 
solidified through free radical polymerization that is initi¬ 
ated by adding 2, 2'-azo-bis-isobutyrynitrile (AIBN) and 
heating to 80 °C for 2 h. The particles then are washed 
and redispersed in deionized water for study. 

Homogeneously nucleated TPM droplets form spheres. 
To synthesize dimpled spheres, we heterogeneously nu¬ 
cleate droplet condensation by adding 0.65 pm-diameter 
polystyrene spheres to the aqueous phase. These small 
spheres serve as nucleation sites for TPM condensation 
and remain embedded in the surface of the resulting 
droplet to a depth that depends on their wetting char¬ 
acteristics. They remain in place during polymerization, 
and thus determine the size of the dimple in the final 
particle. After the TPM is polymerized, the polystyrene 
spheres are dissolved by transferring the particles into 
toluene, leaving uniformly sized dimples. For the sam¬ 
ple represented by Fig. 1(b), the dimple accounts for 5 % 
of the equivalent sphere’s volume. The completed par¬ 
ticles are then transferred back into deionized water for 
cleaning and study. 

Dimpled spheres have immediate applications for lock- 
and-key colloidal self-assembly 7 12 and are models for 
colloidal microcapsules 13 , which are widely used in in¬ 
dustrial applications, and tend to buckle into dimpled 
spheres through osmotic stress 13-15 . They are useful for 


assessing the limits of holographic characterization be¬ 
cause their departure from sphericity is well defined. 

We prepared these particles for holographic character¬ 
ization by dispersing them in water at a volume fraction 
of 10 -5 and introduced into a 50pm-thick gap between 
a glass microscope slide and a cover slip. The sample is 
sealed and mounted on the stage of a custom-built holo¬ 
graphic microscope 1,16 , which illuminates it with a colli¬ 
mated laser beam at a vacuum wavelength of A = 447 nm 
(Coherent Cube). An illuminated particle scatters some 
of the laser light to the focal plane of a microscope ob¬ 
jective lens (Nikon Plan Apo, 100 x, numerical aperture 
1.4, oil immersion) where it interferes with the unscat¬ 
tered portion of the beam. The interference pattern 
is relayed by the objective lens and a tube lens to a 
video camera (NEC TI-324A), which records its inten¬ 
sity at 29.97 frames/s with a calibrated magnification of 
135 nm/pixel. The camera’s 0.1 ms exposure time is short 
enough to avoid blurring of the interference pattern due 
to particle motion 17,18 . The illumination is linearly po¬ 
larized with its axis of polarization aligned to within 1° 
with the x axis of the camera. 

Each video frame is a hologram of the particles 
in the 86 pm x 65 pm field of view. The images in 
Fig. 1(c) and 1(d) are typical holograms of a sphere 
and a dimpled sphere, respectively. They each sub¬ 
tend 150 x 150 pixel and are cropped from the camera’s 
640 x 480 pixel field of view. Such holograms then can 
be analyzed 1,2,17,19 with predictions of the Lorenz-Mie 
theory of light scattering 20,21 to estimate the particle’s 
radius, a p and refractive index, n p . Specifically, we pro¬ 
cess a recorded image /(r) by subtracting off the camera’s 
dark count, Id( r) and normalizing by a background im¬ 
age, Jo(r), that is recorded with no particles in the field 
of view: 


I( r) - I d ( r) 
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Assuming that gradients in the amplitude and phase of 
the illumination are small over the scale of the particle, 
the normalized hologram of a particle located at r p rela¬ 
tive to the center of the microscope’s focal plane may be 
modeled as 1,4 


b(r) = \x + e tkZp t s (k(r - r p ))\ 2 , (2) 

where k = 27rn m /A is the wave number of the light in 
a medium of refractive index n m , and where f s (/cr) de¬ 
scribes how the particle scatters x-polarized light. 

If the particle may be modeled as an ideal isotropic 
sphere, then f s (/cr) is the Lorenz-Mie scattering 
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FIG. 1. (a) Scanning electron micrograph of a colloidal TPM sphere and (b) a dimpled sphere. Scale bars represent 500 nm. 
(c) and (d) Corresponding holograms for particles from the samples show in (a) and (b). (e) 5000 measurements of sphere 
radius and refractive index for a single TPM sphere held in an optical tweezer. Each point represents a single measurement, 
and is colored according to the relative density of measurements, P(a p ,n p ). (f) equivalent result for a dimpled sphere, (g) 
Distribution of size and refractive index for 5000 TPM spheres and (h) 5000 dimpled spheres. 


function 20,21 , which is parameterized by the particle’s ra¬ 
dius and refractive index. Fitting Eq. (2) to a measured 
hologram therefore yields the particle’s three-dimensional 
position r p , its radius a p and its refractive index n p . 
Previous studies on model colloidal spheres have con¬ 
firmed that these fits converge reliably for micrometer- 
scale spheres, and yield the radius with a precision bet¬ 
ter than 5 nanometers 1,3,4 and the refractive index to 
within 3 parts per thousand 3,5 . The data in Fig. 1(e) 
show results obtained for a typical TPM sphere local¬ 
ized in an optical trap that was projected through the 
microscope’s objective lens using the holographic optical 
trapping technique 22 . The spread in values is comparable 
to the numerically estimated uncertainty in the individ¬ 
ual fits, a p = 0.790 =b 0.003 pm and n p = 1.445 =b 0.001, 
suggesting both that the imaging model is appropriate, 
and also that the signal-to-noise ratio estimated by the 
median-absolute-deviation (MAD) metric is reasonable. 
This sphere’s holographically measured radius is con¬ 
sistent with the mean value, 0.76 ±0.06 pm, obtained 
through SEM observations on the same batch of spheres. 
As expected, the polymerized particles’ refractive index 
is larger than that of monomeric TPM oil, 1.431 at the 
imaging wavelength. 

The scattering function for an aspherical object, such 
as a dimpled sphere, depends on the object’s detailed 
shape and orientation 21 . Analytical results are available 
for just a few special cases. For more general cases, nu¬ 
merical methods are required, such as the discrete dipole 
approximation (DDA) 23 . Even with highly optimized 


implementations 24 , however, such approaches are com¬ 
putationally intensive 25 . If an object’s departure from 
sphericity is small enough, and if the influence of the 
non-ideality on the recorded hologram is sufficiently well 
localized within the recorded image, Lorenz-Mie analysis 
may still yield useful results for the particle’s size and 
refractive index without incurring this cost. 

Figure 1(f) shows results obtained by fitting the ideal- 
case model to holograms of an optically trapped dimpled 
sphere. The radius estimated by straightforward Lorenz- 
Mie analysis of 60 such holograms is 0.821 =b 0.006 pm. 
The corresponding estimate for the refractive index, 
1.448 d= 0.001 is remarkably similar to the value obtained 
for the ideal spheres. In both cases, the distribution 
of values extracted from nonlinear least-squares fits to 
Eq. (2) is consistent with single-fit error estimates. This 
agreement suggests that the ideal model can yield quan¬ 
titative results for the characteristics of dimpled spheres, 
without incurring the costs of more realistic modeling. 

SEM analysis suggests that the sample-averaged ra¬ 
dius of the dimpled spheres is 0.75 ±0.05 pm, which is 
significantly smaller than the result obtained holograph¬ 
ically. Similar discrepancies have been noted in previous 
studies 26 , and reasonably may be explained by changes 
induced by preparing the spheres for SEM observation. 

The extent to which a dimpled sphere’s hologram may 
be described with an ideal sphere’s scattering function 
depends on the dimple’s orientation. Optically trap¬ 
ping a dimpled sphere constrains its rotations as well 
as its translations, fixing the dimple’s axis in the trans- 
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verse plane. When released from its trap, the dimpled 
sphere rotates freely in three dimensions, with conse¬ 
quences for holographic characterization. Figure 1(g) 
shows the distribution of characteristics obtained over 
the course of 10 min for a freely diffusing sphere. The 
mean radius and refractive index obtained for this par¬ 
ticle are 0.80 ± 0.01 pm and 1.440 ± 0.003, respectively. 
Additional uncertainty in the particle’s radius and re¬ 
fractive index reflects uncorrected interference artifacts 
in the recorded hologram due to defects in illumination 4 . 
The corresponding distribution for the freely diffusing 
dimpled sphere is broadened still further, and displays 
a strong anticorrelation between radius and refractive 
index. The peak of this distribution remains at the 
trapped-particle values, a p = 0.80 ± 0.01 pm and n p = 
1.447 ±0.005, presumably because the randomly ori¬ 
ented particle is more likely to have its dimple transverse 
to the optical axis than facing it. 

Having the dimple pointing sideways is beneficial for 
holographic characterization. In this orientation, the 
dimple’s contribution to the scattering pattern is asym¬ 
metric, and thus minimally influences the fit to the 
largely symmetric model. It plays a role comparable to 
uncorrected background artifacts, reducing the precision, 
but not seriously affecting the mean values. When the 
dimple is directed along the axis, however, distortions to 
the scattering pattern are symmetric about the axis and 
thus affect the fits more strongly. 

The success of an idealized model for describing light 
scattering by a dimpled sphere may be explained at least 
heuristically by treating the dimple as a volume of the 
sphere whose scattering is phase-shifted by 180°. The 
scattering amplitudes for the sphere as a whole and for 
the dimple scale roughly as their respective volumes. 
If the dimple’s volume is only a small fraction of the 
sphere’s, and if, furthermore, the dimple’s contribution 
is asymmetric, the perturbation should be negligible. 

We applied the same technique to characterizing the 
distribution of properties in monodisperse samples of 
TPM spheres and TPM dimpled spheres. For these mea¬ 
surements, dispersions of particles at a volume fraction 
of 10 -4 were streamed down a channel by a pressure 
driven flow with a peak speed of 100 pm s -1 . This is fast 
enough to acquire holograms of 5000 spheres in 10 min, 
but not so fast as to incur artifacts due to motion-induced 
blurring 17,18 . Results plotted in Fig. 2 for TPM spheres 
reveal a reasonably symmetric distribution of particle 
sizes and refractive indexes peaked at 0.82 ± 0.02 pm 
and 1.446 ± 0.006, respectively. These values are consis¬ 
tent both with the single-sphere measurements reported 
in Fig. 1. The corresponding distributions for dimpled 
spheres plotted in Fig. 2 become increasingly broad and 
asymmetric as the relative size of the dimple increases. 
Particles with 5 % dimple volumes yield a mean refractive 
index, 1.444 ±0.006, consistent with the ideal spheres’ 
value. Increasing the dimple volume to 10% leads to 
substantial deviations, with a mean refractive index of 
1.42 ± 0.03. Some variability can be attributed to the 



FIG. 2. Level sets of the holographically measured distribu¬ 
tions of particle properties for TPM spheres (red), dimpled 
spheres with small dimples (blue) and large dimples (purple). 
The presence of a small dimple has no significant influence 
on holographic characterization results. Larger dimples cause 
systematic errors. SEM images show typical representatives 
of each sample, with scale bars denoting 500 nm. 


dimpled spheres’ random orientation in the channel. 

To assess the extent of the distortions that can be han¬ 
dled with the idealized model for holographic character¬ 
ization, we apply the discrete-dipole approximation 23-25 
to compute holograms of dimpled spheres, and then ana¬ 
lyze the resulting synthetic data with the same software 
used to analyze experimental data. The discrete-dipole 
approximation treats an object as a three-dimensional ar¬ 
rangement of independent microscopic dipole scatterers, 
each of which is illuminated by the incident beam and 
also by the first-order scattering of its neighbors 27 The 
superposition of scattered waves yields an estimate for 
f s (fcr), which is used to synthesize a hologram. Dim¬ 
pled spheres are modeled as the superposition of two 
spheres separated by a center-to-center distance d, one 
of radius a p and refractive index n p: and the other of 
radius and the refractive index of the medium. This 
is shown schematically in Fig. 3(a). Setting = 0 or 
d > a p + dd yields a perfect sphere. Setting d < a p — ad 
with a p > ad yields a sphere with a spherical inclusion. 
Dimpled spheres result when a p — ad < d < a p ±a^. Lines 
in Fig. 3(b) demarcate these regions. 

Analyzing DDA-generated holograms of perfect 
spheres yields excellent agreement with input parame¬ 
ters for sphere radii up to a p < 0.5 pm. Scattering by 
larger spheres requires proper treatment of higher-order 
scattering, and is not supported by the ADDA imple¬ 
mentation of the DDA algorithm that we adopted 24,25 . 
Limiting the analysis to parameters within the spherical 
particle’s domain of applicability, we assessed discrepan¬ 
cies between input parameters and values obtained by 
fitting the resulting holograms with the scattering func¬ 
tion for ideal spheres. 

Figure 3 summarizes the performance of Lorenz-Mie 
analysis for characterizing dimpled spheres. Values in 
Fig. 3(b) represent the factor by which the error in a 








4 



a d [irni] 


FIG. 3. Holographic characterization of dimpled TPM 
spheres with radius a p = 0.5 pm as a function of dimple ra¬ 
dius aa and inset distance, d. (a) Schematic representation of 
sphere geometry, (b) Relative error as a function of dimple 
radius and offset. Colors represent the larger of the errors 
in particle radius and refractive index, relative to their re¬ 
spective estimated uncertainties. Region I represents ideal 
spheres. Region II: dimpled spheres. Region III: spheres with 
spherical inclusions. Region IV: inaccessible. 

dimpled particle’s measured characteristics in increased 
relative to the error for an ideal sphere. As anticipated, 
errors in estimates for a particle’s radius and refractive 
index are smaller than 1 % for dimples that take up less 
than 5 % of the sphere’s volume. Deviations become 
larger for spheres distorted by larger dimples, particu¬ 
larly if those dimples are aligned with the optical axis. 
This leaves a substantial domain of applicability within 
which computationally efficient implementations of holo¬ 
gram analysis can be used to measure the properties of 
this class of imperfect spheres. 
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